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1. Introduction: Flavour Physics vs lattice systematics 

During the last five years, a new generation of experimental facilities dominated by B factories 
have brought Flavour Physics to the era of precision studies. Uncertainties on decay and mixing 
amplitudes of kaons, D- and B-mesons have plummeted, setting more stringent constraints than 
ever on the accuracy required for theoretical estimates aimed at determining Standard Model (SM) 
parameters or probing new physics. Indeed, although a determination of CKM matrix elements 
essentially free from hadronic uncertainties starts to be possible [[[]], knowing all the relevant ma- 
trix elements to the required precision is still essential in order to check the consistency of SM 
predictions and set bounds on effects beyond the SM. 

The requirement of few percent accuracy demands, in particular, a fully first-principles ap- 
proach to the long-distance regime of strong interactions in which all the systematics is consciously 
brought under control. Mandatory features include: 

• Dynamical simulations with at least 2{light) + 1 flavours of sea quarks. 

• Good conceptual control over the regularisation. 

• Good control of all the symmetries (especially: flavour symmetries). 

• Fully non-perturbative renormalisation of all the composite operators involved. 

• Elimination of cutoff dependences. 



Wilson and chirally symmetric fermions are preferred on conceptual grounds. Wilson fermions 
have disadvantages from the point of view of flavour symmetries and treatment of ultraviolet cutoff 
dependencies, but are already able to approach the light quark regime in dynamical simulations 
[§, ||, ^| |5j q|, although recent progress with dynamical chiral fermions has proved equally im- 
pressive (see |0] and references therein). Twisted mass QCD (tmQCD) is a variant of the Wilson 
regularisation that potentially allows for a better control of chiral symmetry breaking and cutoff 
effects, which turns particularly advantageous in the computation of weak matrix elements. In this 
context, it may offer a convenient compromise between an adequate control over chiral symme- 
try and numerical affordability. Precision quenched computations, that properly deal with all the 
systematics apart from dynamical quark effects, are an essential step in order to understand these 
issues and prepare the terrain for dynamical studies. Not least importantly, they offer sound argu- 
ments for the choice of regularisation both for sea quarks and for the valence sector of mixed action 
approaches. 

This paper deals mainly with quenched numerical results for weak matrix elements obtained 
from tmQCD. In Section 2 the ALPHA Collaboration computation of Bk in tmQCD [||] is dis- 
cussed. Section 3 briefly reviews the extension of the strategy to deal with neutral B-meson mixing 
amplitudes, and reports on the project status. Section 4 deals with tmQCD proposals to study 
K — > 7171 amplitudes. Finally, in Section 5 some final remarks are made. 

Much of the work presented |9|, |10[ |ll|] is part of the ALPHA Collaboration research pro- 
gramme. General reviews of progress in kaon and B-physics on the lattice have been provided at 
this conference by W. Lee and T. Onogi [ 12, |13| ]. 
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2. Bk in quenched (tm)QCD 

2.1 Lattice QCD and indirect CP violation in kaon decays 

Indirect CP violation in K — > %% decays is measured by the parameter Sk, defined in terms of 
kaon decay amplitudes as 

_ T(K L ^(nn) I=0 ) 

where / is the total isospin of the two-pion state. Experiment yields [|lj] 

S K = [2.232(7) x lfT 3 ] , e = (43.5 ± 0.7)° . (2.2) 

At leading order in an Operator Product Expansion (OPE) treatment of electroweak interac- 
tions, the Standard Model (SM) prediction for \ek\ can be written as [ ]l5| ] 

\e K \ ~ C e B K lm{V t * d V ts } {Re{V* d V cs } [t 7i 5 (x c ) - mS {x c , Xt )} - Re{V t * d V ts }ji 2 S fa)} . (2.3) 

Here C e = G\F^M K Ml^ / \d\fl% 2 AM K ), Sofa) and S (x c ,x t ) fa = mf/M^) parameterise the Wil- 
son coefficients of the OPE, 1712,3 are short-distance QCD corrections to the latter (known to NLO), 
and 

. (K°\d AS = 2 \K°) 

B K = X ' ' 1 , (2.4) 

where d AS = 2 is the effective four-quark interaction operator 

AS = 2 = (s-^d)(s^d) , (2.5) 

= y^(\ — y 5 ) ? and the hat denotes renormalisation group invariant (RGI) matrix elements. The 
dimensionless parameter Bk thus provides the long-distance, non-perturbative QCD contribution, 
and largely dominates the uncertainty on the SM value for \ek\. In the standard Unitarity Triangle 
(UT) analysis of CP violation in the SM, the value of |eg-| provides a hyperbola in the (p,fj) plane. 
After the recent generation of experimental results from B-factories, this is one of the least precise 
UT constraints. Improving the accuracy of Bk is hence essential in order to derive stringent bounds 
on the amount of non-SM CP violation in kaon decay. 

Besides quenching, which is an uncontrolled source of systematic error, the most important 
source of uncertainty in lattice QCD computations of Bk with Wilson fermions arises from operator 
renormalisation. In standard notation, the operator o AS=2 is customarily split into parity-even and 
parity-odd parts as 

qAS-2 _ Q yy+AA _ o VA+AV . (2.6) 

Since parity is a QCD symmetry, the only contribution to the K Q -K Q matrix element comes from 
Ovv+aa- In regularisations which respect chiral symmetry, the latter operator is multiplicatively 
renormalisable. If chiral symmetry is not preserved, C?vv+aa mixes with four other dimension-6 
operators [|l6|, |7|, [18[ [l9|, with positive parity: 



(0 r )vv+aa(Ai) =Z V v+aaC?o,«M) 0vv+AA(go) + £A ; (g o )0,-( 

1=1 



4 



(2.7) 
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The operators 0,-(go) belong to different chiral representations than C?vv+aa- The mixing coeffi- 
cients A,- (go) are finite functions of the bare coupling, while the renormalisation constant Zw+aa 
diverges logarithmically in afX. 

Two proposals have attempted to eliminate operator mixing. They are both based on the ob- 



servation [18, 20] that, even in the absence of chiral symmetry, the operator C?va+av is protected 



from finite operator mixing by discrete symmetries, and thus it renormalises multiplicatively, viz. 

(Or)va+av(a0 = Z V A+Avteo,fl'A t )OvA+Av(go) • (2.8) 

The first proposal [^l|] consists in obtaining the physical K° — K° matrix element of Ovv+aa from 
a correlation function of the renormalised operator C?va+av> related to it through axial Ward iden- 



tities. The method has been put to test in ref. Q22[], with the result that the Bk estimate turned out 



to be compatible with the result of computations that involve operator subtractions. Unfortunately, 
the correlation function of C?va+av is a four-point function, while the matrix element of C?vv+aa 



can be extracted from a three-point function. Thus, the conclusion of [22] is that this method is 
successful in eliminating an important source of systematic errors (operator subtraction) at the cost 
of increased statistical fluctuations. 



In the work under consideration here, the second proposal [)23[], based on twisted mass QCD, 
is implemented. In tmQCD the breaking pattern of flavour symmetries is controlled by the value 
of the twist angle; in particular, the latter can be tuned so as to preserve part of the axial subgroup, 
at the price of breaking vector symmetries, as well as parity. It is thus possible to set up regularisa- 
tions in which the renormalisation of composite operators is greatly simplified. The relevant case 
for us is the renormalised (K°\0\]\j + \\\K Q ) matrix element, which via the tmQCD formalism can 
be extracted from a three-point conelation function of the operator C?va+av- As the tmQCD action 
differs from the standard Wilson fermion action by a soft term, the renormalisation properties of 
composite operators in mass independent renormalisation schemes are not modified. In particular, 
Ova+av remains multiplicatively renormalisable, with the same renormalisation constant and run- 
ning as with Wilson fermions. Thus finite subtractions are avoided in the tmQCD determination of 
B K . 

An obvious alternative to avoid renormalisation problems consists in using regularisations 
with exact chiral symmetry. However, the computational costs involved make it difficult to perform 
continuum limit extrapolations and study finite volume effects. 1 In the case of staggered fermions, 
apart from the operator mixing (the details of which depend on the specific setup), some additional 
problems are present — large scaling violations unless high levels of 0(a 2 ) improvement are im- 
plemented, uncertainties related to the choice of interpolating operator, as well as other difficulties 
related to the breaking of flavour symmetries and the presence of unphysical flavours. 2 Wilson 
fermions therefore offer, a priori, a good compromise between good control of the field-theoretical 
aspects of the problem and affordable computational costs. 



For a state-of-the-art determination of AS = 2 matrix elements, see | 
2 See HO] for an updated discussion of staggered quark results. 
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2.2 tmQCD setup 

We will employ two different fermion actions, namely 

S { * ,2) =a 4 £ [xjf(x)(D w>m +mi + ir 5 T: 3 Hi)Y(x) + s(x)(D w , sw + m s )s(x)] , (2.9) 
x 

4 w/4) =a 4 £ [u{x)(D w ^ + m u )u{x) + ^(x)(D w , sw + m, + //sT 3 ^)^*)] . (2.10) 

The labels on the two actions refer to the values that will eventually be set for the twist angle. 
In Eq. (13) y = (u,d) T , while in Eq. ( 2.1C ) y = (s,d) T . In both cases, the matrix T 3 acts on 



flavour space, and D w sw is the Wilson-Dirac operator with a Sheikholeslami-Wohlert term; c sw is 



tuned to its non-perturbative value [25]. In Eq. (2.1C) it has been assumed a priori that the s and d 



quarks have degenerate physical masses; while this is not necessary as long as this action is used 
in quenched QCD, all the computations carried out with it are performed in that limit. The action 



in Eq. (2.9), on the other hand, is perfectly well suited for an unquenched computation, and it has 



been used to explore the effect of having non-degenerate s and d quark masses (see below). 



The properties of tmQCD have been extensively discussed in several publications (see [23 



26, 27] and references therein). Here we just remind some basic facts. The physical renormalised 



masses of twisted quarks and the twist angle a are given by 



M R = V m R + ' U R' (2 - 11} 

tan a = —, (2.12) 
m R 

where tmr (resp. ;Ur) are the renormalised standard (twisted) quark masses. In order to tune the 
twist angle to some prescribed value up to 0(a 2 ) corrections, we employ the formulae for the 
construction of 0(a) improved renormalised masses 

m R j =Z m [m q j(\+b m am q j) + b m aixf], (2.13) 
HR,l=Z^i(l+b^am q; i), (2.14) 

where m q j = |(l/fc — l/fQ^isthe subtracted bare standard quark mass. 

In the case of the 7t/2 regularisation, in order to have a = n/2 it is enough to set tmr / to zero, 
which is achieved by setting 

flm,,/ = -b m (ajXi) 2 . (2.15) 

The 7r/4 case is somewhat less trivial. Setting a = n/4 requires jIIr,/ = m^;, which via Eqs. 
( |2.13| , [2.14[ ) translates into 



am q ^i = ^-ajXi { 1 + [— (bfj, - b m ) - ZZ A b m ] a\l\ } (2. 16) 

with Z = Z m j (Z^Za). For a given choice of a/i/, K is tuned so that am q j satisfies one of the two 
above relations, taking the values of fc c and all the renormalisation constants and improvement co- 
efficients involved as input. The precision to which the latter are known poses an implicit constraint 
on the accuracy of the tuning of the twist angle. 
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Contact with QCD is made via the change of fermion variables 

y ^ \f/ = R(a)\if, iff -> iff' = \jfR(a) , (2.17) 

where R(a) = expiree! 3 } and y is the twisted quark doublet. This axial rotation induces a 
mapping between composite operators in tmQCD and QCD, which is realised at the level of renor- 
malised correlation functions (or, alternatively, renormalised matrix elements). The relation we are 
most interested in is 

{K°\ (Ofi)w+AA |^°)qCD = — i(K | (Or)vA+AV |^°)tmQCD , (2.18) 

which holds in the continuum limit for the two versions of tmQCD under consideration. From this 
identity, Bk can be extracted from a K°-K° matrix element of the multiplicatively renormalisable 
operator C?va+av- 

It is important to stress that none of the above setups leads to a computation of Bk that involves 
fully twisted quarks only. Hence, the automatic 0(a) improvement argument of Frezzotti and Rossi 
[ ^8[ ] does not apply, and in order to have full 0(a) improvement of the matrix element it would be 
necessary to subtract a number of dimension-seven counterterms from the four-fermion operator. 
Such a procedure is highly impractical, and has not been pursued. Hence, leading cutoff effects in 
Bk are expected to be linear in a. 

2.3 Renormalisation 

The non-perturbative renormalisation of the operator C?va+av has been addressed in ^ [lO|] us- 
ing standard Schrodinger Functional (SF) techniques (see e.g. J29|]). After having defined suitable 
SF intermediate renormalisation schemes, a recursive step-scaling procedure allows to compute to 
high accuracy the renormalisation group (RG) running of the operator in quenched QCD in the con- 
tinuum limit from a low-energy reference scale (2L max ) _1 close to Aqcd to scales of 0(100 GeV), 
where reliable contact with perturbation theory can be made. Together with the renormalisation 
constants at (2L max ) _1 , this provides RGI renormalisation factors free from any uncontrolled sys- 
tematic uncertainty. 

In the particular case of C?va+av> nine different SF schemes were defined and found to provide 
consistent results for RGI renormalisation factors. The continuum limit of the RG running was 
controlled by performing independent simulations with two different fermion actions (plain and 
0(a) improved Wilson fermions). The quality of the result is illustrated by Figure [l]. This approach 
is currently being pursued in order to extend the non-perturbative renormalisation of Ova+av to 



Nf = 2 QCD [J30J]. It is also worth mentioning that the scope of [g, |10j goes well beyond the case 



of the AS = 2 effective Hamiltonian. For instance, [ pl[ ] made use of the results in [g] to address the 
renormalisation of the AS = 1 effective Hamiltonian with an active charm quark constructed with 
overlap fermions. To that purpose, the logarithmically divergent renormalisation constants required 
have been computed through a matching of non-perturbatively renormalised RGI tmQCD matrix 
elements to bare overlap matrix elements at a reference mass mps ^ rriK ■ This procedure is similar 
to the one employed in [[52|] for the renormalisation of the quark condensate. 
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Figure 1: RG running (9va+av in quenched QCD in the SF scheme. Non-perturbative values (circles) are 
compared to perturbation theory predictions. 



2.4 Simulations for bare matrix elements and systematics 

The bare value of Bk can be extracted from the ratio of SF correlation functions 

3 — (/VA+Av(*()) 



R{xq) 



16 [Z A (/ A (^ ) +c A ad f P (x Q )) - iZyf y (xo)][Z A (fl(xo) + c A adof£(xo)) - /Z v /^(x )] ' 

(2.19) 



where /va+av is the correlation function of C?va+av with two pseudoscalar SF boundary sources 
and f x (X = Ao,P,Vq) is a two-point function of the bilinear operator X with a pseudoscalar SF 
boundary source. Precise definitions can be found in [0]. The combination of currents in the 
denominator corresponds to the physical, O(a) improved, renormalised axial current via the chiral 
rotation in Eq. ( 2.17| ). 

Quenched simulations have been performed at j3 = 6.0,6.1,6.2,6.3 for the n/2 regularisation 
and j6 = 6.0,6.1,6.2,6.3,6.45 for n/4. The physical masses of the s and d quarks have always been 
kept degenerate, save for a subset of simulations meant to probe SU(3) flavour breaking effects (see 
below). In the n/2 case, since the s quark is untwisted, it is impossible to reach pseudoscalar masses 
in the region of m^, due to the presence of exceptional configurations. Therefore, pseudoscalar 
masses larger than tuk are simulated, and the results are hence extrapolated to the physical kaon 
mass. In the n/4 case, on the contrary, values at tuk can be obtained by interpolation. 3 The 
necessary renormalisation constants and improvement coefficients, as well as the values of jq, 
have been gathered from the literature. 4 The scale is always fixed via the ratio r^/a as given by 
Q, with r = 0.5 fm. 

The results have been subjected to a number of checks, meant to assess various systematic 
uncertainties: 



• Finite volume effects. In the 7i/2 case, simulations at /3 = 6.0 have been performed for 
physical lattice sizes around L 1.5 fm and L ps 2 fm, at the lowest pseudoscalar meson 

3 The ft = 6.45 data, however, have been obtained at larger masses and then extrapolated to m.K ■ 
4 See Appendix A of [M] and references therein. 
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mass available. Consistent values for all the relevant observables are obtained within errors. 
Simulations were hence carried out on lattices with L ~ 1.5 fm. In the n/4, a similar study 
was performed at j8 = 6.0,6.2 and masses around m^. In this case the conclusion is that 
lattice sizes L ~ 2 fm are needed to avoid finite volume effects. 

SU(3) breaking effects. Physical SU(3) breaking effects on B& have been studied in the n/2 
case at /3 = 6.0. To this purpose, simulations were performed for three values of the ratio £ = 
(M s — Md) / (M s + Md) (where M f is the physical quark mass), namely £ = 0.00,0.16,0.41, 
at fixed romps = 1.78. No effect was observed on Bk within uncertainties, hinting at small 
SU(3) breaking effects. It has to be stressed, however, that the simulated pseudoscalar mass 
is relatively high. 

Spurious SU(3) breaking. The 0(a 2 ) breaking of vector flavour symmetries induced by the 
presence of the twisted mass term has received considerable attention in the literature (see 



[26]). In order to check its effect, one can compare the pseudoscalar meson mass obtained 
for various flavour combinations of twisted and untwisted quarks in the n/2 case. In the 
M s = Md limit, the resulting states can be interpreted as belonging to a multiplet of pseudo- 
Goldstone bosons; hence, deviations from unity in the ratios (mps/rap S ) 2 of squared pseu- 
doscalar masses in different flavour channels quantify the 0(a 2 ) vector symmetry breaking. 
The values of these ratios show that the effect is never beyond the few percent level, and con- 
verges to zero in the continuum limit. The splitting tends to grow mildly as the quark mass is 
decreased. These findings contrast, but are by no means incompatible, with the observation 



of much larger amounts of 0(a) flavour breaking at lighter quark masses [34]. 



After the publication of a more detailed analysis of the accuracy of the tuning of quark 
masses and twist angles was performed. The quality of the tuning was found to be satisfactory in all 
cases save for the simulations at j8 = 6.1, mainly in the n/4 case. This is signalled e.g. by relatively 
large differences between the value of the target twist angle, set to n/2 or n/4 when tuning the 



quark masses via Eqs. ( |2.13U2.14| ), and the value obtained by computing the ratio /iR,//mR / with 
the PCAC quark mass instead of the subtracted quark mass. 

The reason for this behaviour has been traced back to the value of fc c taken as input from the 
literature. Indeed, for an accurate determination of it is crucial to follow a constant physics 
condition in the approach to the continuum limit, which fixes the 0(a 2 ) ambiguities coming from 
this source. Instead, the value = 0.135496 quoted in [[35|] comes from an interpolation of data 
obtained from a constant physics condition at other values of j3. While the effect of relaxing the 



constant physics requirement was found to be negligible for the data of [35], its impact on the tuning 
of twist angles is large. The j8 = 6.1 critical point has been hence determined afresh, obtaining 
Kc = 0.135665(11), and j3 = 6.1 simulations with new mass parameters have been performed. Full 



details will be provided in a forthcoming publication [|36j] . 
2.5 Continuum limit 

As noted above, taking the continuum limit for B% involves a linear extrapolation in a. At this 
stage, having results from two different regularisations, which can be combined in a fit constrained 
to a common continuum limit, is essential for a proper control of the extrapolation. 
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Figure 2: Uncertainties on fi/f related to the O(a) improvement of bilinears. The left panel displays 7T/2 
data, while the right panel shows n/4 data. 



It turned out that one of the most relevant sources of cutoff effects is related to the construction 



of the O(a) improved bilinears in the denominator of Eq. ( 2.19 ). For instance, using either the 
values for Z A ,Z v ,c A determined by the ALPHA Collaboration or those obtained by the LANL 
group [37] results in sizeable effects on Bk at the lowest values of j3 available (see Figure^]). This 
signals the presence of large 0(a 2 ) ambiguities in Bk far from the continuum limit. Combined 
linear+quadratic extrapolation of the data proved to be unstable. Thus the values of j8 for which 
the difference between ALPHA and LANL constructions of O(a) improved bilinears results to Bk 
discrepancies beyond one sigma were conservatively discarded in the linear fits to the continuum 
limit. This means that results at j8 = 6.0 and j8 = 6.1 had to be left out. The resulting extrapolation 
is illustrated by the left panel of Figure 0. The final results are: 



Bx = 0.735(71) . 



fi^ s (2GeV) =0.534(52). 



(2.20) 
(2.21) 



When comparing with the result quoted in [Q], it has to be taken into account that /3 = 6.1 data 
have been revised, for the reasons explained above. 

The value for Bk in Eq. ( 2.20 ) is shown in the right panel of Figure || alongside other repre- 
sentative results in quenched QCD found in the literature. As discussed in Appendix E of [g], the 
difference with other computations with Wilson fermions is mainly due to the method employed to 
determine Bk'. instead of using a ratio similar to the one in Eq. ( 2.19D , the authors of [^] extract 
Bk from a fit of the mass dependence of a different ratio of correlation functions, inspired by Chiral 



Perturbation Theory. It has to be stressed that the computation of [ 22 ] does not have direct access 
to the physical kaon mass region. 

The result of Eq. ( 2.20Q is the only existing quenched result in the literature which has si- 
multaneously eliminated any systematic uncertainty related to renormalisation (both at a reference 
scale and from the point of view of RG running), ultraviolet cutoff dependences, and finite volume 
effects (within the available accuracy). On the other hand, the control of the mass dependence of 
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Figure 3: Left: Continuum limit extrapolation of Bk- Right: Comparison with other quenched results. 



B K with Wilson fermions is still not as accurate as with e.g. Neuberger or domain wall fermions. 



Overall, it seems fair to claim that Eq. ( |2.20| ) is a benchmark result for Bk in quenched QCD. 



3. A strategy to compute B B 

Long-distance QCD contributions to indirect CP violation in the B-meson sector of the SM are 
encoded in the bag parameters 



0|/)AB=2|gO 



Be, 



(B u |0; 



where the relevant effective four-quark interactions have the form 

(bp)(bp), £ = d,s. 



(3.1) 



(3.2) 



These B-parameters appear, together with the corresponding meson decay constants, e.g. in the 
expressions for neutral B-meson mass differences AM( [ |L5| ] 



AM d = 0.05 ps~ 



AM S = 17.2 ps 



-l 




(3.3) 
(3.4) 



where tjb encodes short-distance QCD effects. The experimental values for these quantities are 
AM d = 0.507 ± 0.005 ps" 1 pi] and AM,. = 17.77 ± O.lO(stat) ± 0.07(sys) ps" 1 @. The recent 



measurement of AM S by the CDF Collaboration has set very stringent constraints on the required 
precision of theoretical determinations of Bs s Fg , which are now at the same level as those on 
BB d Fg . In addition to "standard" systematic uncertainties such as dynamical light quark effects, 
matrix elements involving heavy quarks are particularly sensitive to improvements coming from a 
systematic, conceptually controlled treatment of heavy quark effects within lattice QCD. 
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A strategy for a precise lattice QCD computation of Bs n which in principle would keep all 



systematic uncertainties under control, has been put forward in JllJ]. The b quark is treated at 
leading order in Heavy Quark Effective Theory (i.e. in the static approximation), although I /mi, 
corrections from the heavy quark expansion can be eventually included following the spirit of 



[p9[]. In the static approximation the relevant physical amplitude is a linear combination of matrix 
elements of two static-light four-fermion operators, viz. 

Of =2 = (h^l)(h^£), l = d,s, (3.5) 
Of= 2 = (h(l- K)£)(h(l-K)£), £ = d,s. (3.6) 

The renormalisation of generic static-light four-fermion operators with Wilson light fermions has 
been analysed in detail in [O]. An important conclusion of this study is that, similar to the case of 
fully relativistic operators, parity-even operators mix between them due to the breaking of chiral 
symmetry. On the other hand, in the parity-odd sector it is possible to find a complete basis of op- 
erators that renormalise multiplicatively. This opens the door to a generalisation to this context of 
the tmQCD strategy pursued for Bk ■ In particular, it is possible to extract Bb, from matrix elements 



of the VA + AV and SP + PS parts of the operators in Eqs. ( |3.5| , p.6| ) if the light quark flavour I is 
twisted at a = n/2. This avoids any need of dealing with complicated operator renormalisation pat- 
terns, and eliminates any constraint on quenched computations due to exceptional configurations. 
Furthermore, it is possible to extend the automatic 0(a) improvement arguments of Frezzotti and 
Rossi to show that the matrix elements of interest will display scaling violations at 0(a 2 ) only. 

The numerical implementation of non-perturbative renormalisation for static-light four-fermion 
operators is discussed in detail in [fill]. Preliminary results for the RG running of static-light four- 



fermion operators in quenched QCD are shown in Figure H. Final results will be the object of a 



forthcoming publication [4C]. 



4. tmQCD for K -> %%1 

The computation of non-leptonic kaon decay amplitudes in QCD poses much harder problems 
than those related to AF = 2 processes: 

• Finite volume effects strongly affect the two-pion final state, making the direct extraction 
of the amplitudes from Euclidean correlation functions considerably difficult pT] , |42| ]. It has 
thus become customary to attempt instead the computation of the relevant couplings in a low- 
energy effective description of QCD based on Chiral Perturbation Theory, which in principle 
would allow the computation of the amplitudes at a given order in the chiral expansion p3|]. 
This requires, ideally, access to the chiral regime of QCD. 

• The renormalisation of the AS = 1 effective Hamiltonian arising from the OPE treatment 
of electroweak interactions requires dealing with a complex operator mixing problem. In 
particular, if chiral symmetry is not preserved by the lattice regularisation, as with Wilson 
fermions, mixing with lower dimension operators proceeds via coefficients that diverge with 



an integer power of the cutoff [44]. On the other hand, if the charm quark is kept as an 
active degree of freedom the presence of exact lattice chiral symmetry eliminates all power- 
divergent mixings [45, |46| ]. 
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Figure 4: Preliminary results for the non-perturbative RG running ("M.C." points) of the four independent 
elements of the multiplicatively renormalisable basis for parity- odd static-light four-fermion operators. See 
[ ITU for notational details. The solid curves are perturbative predictions. 



The case for employing regularisations that preserve chiral symmetry in the approach to this 
problem is therefore very strong. Indeed, chiral fermions have been instrumental in a recent com- 
putation, in the quenched approximation, of the leading-order low-energy couplings of the AS = 1 
effective weak Hamiltonian in the GIM limit m c = m u [|47|], the first results of a comprehensive 
programme aimed at understanding the role of the charm quark in the AI = 1/2 enhancement rule 
[46, 48]. On the other hand, it is conceivable that the control over chiral symmetry breaking that 
constitutes one of the main assets of tmQCD can be exploited in order to alleviate the problems 
related to renormalisation. Two different proposals have been actually put forward to that effect. 

In [49], a theory with four quark flavours is considered, without specifying a priori how many 
of them are dynamical. Once the light doublet is twisted at angle tt/2, it is immediate to show that 
the power divergences affecting K — » % matrix elements are at most linear, and that there are no 
finite mixings with other dimension six operators. This is a substantial gain with respect to standard 
Wilson fermions, in which divergences are quadratic and finite mixings are present. If the heavier 
s,c flavours are fully twisted as well (which is straightforward if they are kept quenched), then it is 
possible to eliminate power divergences altogether, simply by employing non-perturbatively 0(a) 
improved fermion action and quark bilinears. 



In Q50|], the authors consider a theory in which a valence sector containing an arbitrary number 
N v of flavours is matched to a theory with Nf dynamical quarks. All the flavours are fully twisted. 
The freedom to fix twist angles arbitrarily for valence quarks, without the need to restrict to non- 
anomalous chiral rotations, is then used to set up a valence sector that allows to extract K — > % 
matrix elements from correlation functions that do not require any power divergent subtraction. 
The authors propose a specific valence sector with N v = 10. In the same paper, a similar technique 
is proposed to obtain a multiplicatively renormalisable Bk\ in this case, N v = 6. A strong advantage 
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of this framework is that only fully twisted quarks are used, and hence automatic 0(a) improvement 
arguments apply. 

These tmQCD proposals are appealing in that they potentially offer many of the advantages of 
exactly symmetric regularisations at a considerably lower computational cost. It has to be stressed, 
however, that the arguments which show that undesired counterterms cancel rely crucially on the 
assumption that a precise tuning of the twist angle has been performed. In the case of power di- 
vergences the issue is particularly sensitive, as systematic uncertainties in the tuning of parameters 
may result in a lack of cancellation of large contributions to correlation functions. It is important to 
notice, too, that the absence of exact chiral symmetry poses an intrinsic lower bound to the quark 
masses that can be simulated safely; in particular, access to the deep chiral regime, as achieved 
in p7|], may be compromised. Finally, the need to separate the AI = 3/2 and AI = 1/2 channels 
requires a good control over the 0(a 2 ) breaking of isospin symmetry inherent to tmQCD. Given 
these caveats, the suitability of tmQCD to deal with K — ► Tin decays is an open problem that may 
only be settled by dedicated numerical studies. 



5. Conclusions 

Twisted mass QCD, together with state-of-the-art techniques for Wilson fermions, allow for 
benchmark quenched computations of weak matrix elements, as shown by Bk- The ideas put 
forward for AS = 2 matrix elements can be extended to other problems, like AB = 2 and K — > Tin 
amplitudes, offering potential for precise computations that do not resort to exact chiral symmetry. 

The dominant source of uncertainty left in the quenched approximation (certainly so for Bk) 
is related to the lack of full 0(a) improvement, which amplifies the error of the continuum limit 
extrapolation. Thus, if Wilson fermions are to be used in the future in the determination of weak 



matrix elements, the use of tmQCD variants that embody automatic 0(a) improvement [5C] may 
prove essential. Two important aspects of the tmQCD approach are critical in the context of weak 
matrix elements: the tuning of parameters, in particular of the twist angle, has to be controlled 
to high precision; and flavour symmetry breaking effects have to be kept at the few percent level. 
The question whether valence tmQCD quarks offer a convenient alternative to chirally symmetric 
fermions for some specific applications remains to be addressed by dedicated simulations. 
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